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The paper deals with the Neumann spectral problem for a singularly perturbed second 
order elliptic operator with bounded lower order terms. The main goal is to provide a 
refined description of the limit behaviour of the principal eigenvalue and eigenfunction. 
Using the logarithmic transformation we reduce the studied problem to additive eigenvalue 
problem for a singularly perturbed Hamilton-Jacobi equation. Then assuming that the 
Aubry set of the Hamiltonian consists of a finite number of points or limit cycles situated 
in the domain or on its boundary, we find the limit of the eigenvalue and formulate 
the selection criterium that allows us to choose a solution of the limit Hamilton-Jacobi 
equation which gives the logarithmic asymptotics of the principal eigenfunction. 
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1 Introduction 


This paper is devoted to the asymptotic analysis of the first eigenpair for singularly perturbed 
spectral problem, depending on the small parameter e > 0, for the elliptic equation 


Ettij (x) 


d'^Ue 

dxidxi 


r / \dUe , , , 

+ = XeU 


( 1 . 1 ) 


in a smooth bounded domain hi C with the boundary condition 


du^ 

du 


= 0 , 


( 1 . 2 ) 


on do,, where ^ denotes derivative with respect to the external normal. 

The bottom of the spectrum of elliptic operators plays a crucial role in many applications. In 
particular, the first eigenvalue and the corresponding eigenfunction of fll.ip - fll.2p . are important 
in understanding the large-time behavior of the underlying non-stationary convection-diffusion 
model with reflecting boundary. Due to the Krein-Rutman theorem the first eigenvalue Ag of 
fll.ip - fll.2p (the eigenvalue with the maximal real part) is simple and real, the corresponding 
eigenfunction Us can be chosen to satisfy Us{x) > 0 in D. 

The goal of this work is to study the asymptotic behavior of and Mg as e —)■ 0. While 
in the case of constant function c{x) in (II.ip the first eigenpair is (trivially) explicitly found, 
the asymptotic behavior of the first eigenpair is quite nontrivial when c{x) is a nonconstant 
function, in particular, the eigenfunction might exhibit an exponential localization. 

Boundary-value problems for singularly perturbed elliptic operators have been actively stud¬ 
ied starting from 1950s. We mention here a pioneering work IZB. where for a wide class of op¬ 
erators (so-called regularly degenerated operators) the asymptotics of solutions were obtained. 

In the works [19], [20|, [6] (see also [5]) the principal eigenvalue of singularly perturbed 
convection-diffusion equations with the Dirichlet boundary condition was investigated by means 
of large deviation techniques for diffusion processes with small diffusion. In [2] the estimates 
for the principal eigenvalue were obtained by comparison arguments and elliptic techniques. 

The case when convection vector field has a finite number of hyperbolic equilibrium points 
and cycles was studied in [8] where methods of dynamical systems are combined with those 
of stochastic differential equations. These results were generalized in [3] to the case when the 
boundary of domain is invariant with respect to convection vector held. Similar problem in the 
presence of zero order term was considered in HD]. 

In [15] the viscosity solutions techniques for singularly perturbed Hamilton-Jacobi equation 
were used in order to study the principal eigenfunction of the adjoint Neumann convection- 
diffusion problem. The logarithmic asymptotics of the eigenfunction was constructed. 

The work [T3] deals with the principal eigenpair of operators with a large zero order term 
on a compact Riemannian manifold. The approach developed in this work is based on large 
deviation and variational techniques. 
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Dirichlet spectral problem for a singularly perturbed operators with rapidly oscillating lo¬ 
cally periodic coefficients was studied in [16] and m- In [16] with the help of viscosity solutions 
method the limit of the principal eigenvalue and the logarithmic asymptotics of the principal 
eigenfunction were found. These asymptotics were improved in ng and HZ] using the blow up 
analysis. 

In the present work when studying problem fll.ip - fll.2p . we make use of the standard vis¬ 
cosity solutions techniques in order to obtain the logarithmic asymptotics of the principal 
eigenfunction. However, the limit Hamilton-Jacobi equation in general is not uniquely solvable 
and does not give information about the limit behaviour of A^. Therefore, we have to consider 
higher order approximations in (ll.ip - (ll.2p . Under rather general assumptions on the structure 
of the Aubry set of the limit Hamiltonian, we hnd the limit of Ag and can choose the solution 
of the limit problem which determines the asymptotics of the principal eigenfunction. Notice 
that we did not succeed to make the blow up analysis work in the case under consideration. In 
this case, for components of the Aubry set located on the boundary, the natural rescaling still 
leads to a singularly perturbed operators . Instead, we study a rehned structure of solutions of 
the limit Hamilton-Jacobi equation in the vicinity of the Aubry set. This allows us to construct 
test functions that satisfy the perturbed equation up to higher order. 

We also would like to remark that, with obvious modihcations, the results of this work as 
well as the developed techniques remain valid for the boundary condition of the form 


du^ 

<9/3 


= 0 , 


where /3 is a U^-smooth vector held on dVL non-tangential at any point of dVL. In particular, 
conormal vector held = aijPj can be considered. 


2 Problem setup and results 

We study problem (ll.ip - (ll.2p under the following assumptions on the operator coefficients and 
the domain: 

(al) n is a bounded domain in N > 2, with a boundary; 

(a2) all the coefficients are U^-functions in H; 

(a3) the matrix [aij) is symmetric and uniformly elliptic. 

Further assumptions on the vector held b will be formulated later on. 

Since > 0 in hi we can represent in the form 
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this results in the following nonlinear PDE 


1 

- ttij [x) ——-h -H (Vhhe, x) + c{x) = Ae in 

oxiOXj e 


or 


d‘^W, 

- Hiyw^.x) +ec{x) = eXs 

OXiOXj 


in 


with the boundary condition 


where 


dWe 

du 


= 0 on dVL^ 


H[p,x) ^ aij{x)piPj - bi{x)pi 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 

(2.4) 


is a function to be referred to as a Hamiltonian. Passing to the limit as £ —1 0 in (12.211 . with the 
help of the standard approach based on the maximum principle, we can show that We converges 
uniformly (up to extracting a subsequence) to a viscosity solution W(x) of the Hamiltion-Jacobi 
equation 

H{VW{x),x)=0 infJ (2.5) 

with the boundary condition 

dW 

= 0 on (2.6) 

ou 

Recall that a function W ^ C (fJ) is called a viscosity solution of equation (12.bh if for every 
test function $ G C'°°(r2) the following holds 


• if hh — $ attains a maximum at a point G hi then W (V$(.^), .^) < 0; 


• iiW — ^ attains a minimum at G hi then hh(Vd)(,^), ,^) > 0 . 

The boundary condition fl2.6p is understood in the following sense, V<h G C'°°(r2) 

• if IT — $ attains a maximum at G dVt then min|jf(V<h(.^), .^), fy(0| — Oj 

• if IT — $ attains a minimum at ^ G dVt then max|iJ(VT(^), ^), fp(0| — 0- 
It is known [7] that every solution of problem fl2.5p - fl2.6p has the representation 

W{x)= ini^[dH{x,y) + W{y)y (2.7) 

where Ah is so-called Aubry set and y) is a distance function. To define Ah and dnix, y) 
consider solutions of Skorohod problem 

' _ 

77(f) G r 2 ,f > 0 

< 77(f) + a{t)p{ri{t)) = v{t) with a{t) > 0 and Q;(f) = 0 when 77(f) ^ dO, (2.8) 

77(0) = X 
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where v G -h^((0, cxo);R^) is a given vector field and a; G is a given initial point, while the 

curve 7] G oo); and the function a G L^((0, oo); M+) are unknowns. Under our 

standing assumptions on (c)f2 G U^) Skorohod problem fl2.8p has a solution, see [7]. 

Consider now the Legendre transform L{v,x) = sup {v -p — H{p,x)) and define the distance 

pgR^ 

function 


dH{x,y) = inf < / L{—v{s),p{s))ds, p solves fl2.8p . p{0) = x, p{t) = y, t>0 


(2.9) 


Next we recall the variational definition of the Aubry set 


xeAH^^^yS>0 77 solves ( 12 . 8 p . pfO) = p(t) = x. f > = 0 . 

( 2 . 10 ) 

In this work we assume that the Aubry set has finite number of connected components 


Ah = \^ Ak and each Ak is either isolated point 

finite 

or closed curve lying entirely either in or on dQ. ( 2 . 11 ) 

Additionally we assume that 


if Ak Q then Ak is either hyperbolic fixed point 

or hyperbolic limit cycle of the ODE x = b{x); (2.12) 


if Ak C dfl then the normal component bu{x) of the field b{x) is strictly positive on Ak 
and Ak is either hyperbolic fixed point or hyperbolic limit cycle of the ODE 
X = b^{x) on dD, where b^-^x) denotes the tangential component of b{x) on dQ. (2.13) 

In order to state the main result of this work we assign to each component Ak of Ah a 
number a{Ak) as follows. If Ak is a fixed point of the ODE x = b{x) and ^ G D, linearizing 
the ODE near ^ to get z = B{ff)z we define a{Ak) by 

a{Ak) = -J2^^ + <0, (2.14) 

9i>0 

where 6i are the real parts of eigenvalues of the matrix B{^). Note that the hyperbolicity of 
the fixed point means that the eigenvalues of B{^) cannot have zero real part. If Ak = {O and 
^ G dQ, consider the ODE x = b^-^x) on dD in a neighborhood of the point Passing to the 
linearized ODE z = Br{^)z in the tangent plane to dQ at the point we denote by 9i the real 
parts of the eigenvalues of i?r (0 and set 

a(A) = -5^^'* + c(0, (2.15) 

h>o 


5 






Consider now the case when {Ak} C is a limit cycle of ODE i. = b{x). Let P > 0 be the 
minimal period of the cycle and let ©j be the absolute values of eigenvalues of the linearized 
Poincare map. (Recall that the limit cycle is said hyperbolic if there are no eigenvalues of 
linearized Poincare map with absolute value equal to 1.) We dehne now cr{Ak) by setting 


(T{Ak) 


^^log0i + ^/ c{^{t))dt, 

^ 0i>i 


(2.16) 


where ^{t) solves ^ ^{t) G Ak- 

Finally, in the case when > 0 on Ak and Ak is a limit cycle of the ODE x 
we set 


a{Ak) = 


0i>i 


c(^(f))df, 


hr{x) on (90, 

(2.17) 


where ^ = fer(0 ^md ^{t) G Ak, P is the minimal period and ©j are the absolute values of the 
eigenvalues of the linearized Poincare map. 

The main result of this work is 


Theorem 1. Let conditions (aA)-(aA) be fulfilled, and assume that the Aubry set Ah satisfies 
\2.11\) . and h2.1A) . Then the first eigenvalue of U.l\) converges as e ^ 0 to 


linr Ap 

e -)-0 


max|(7(.4fe); Ak C Au 


(2.18) 


where a{Ak) is given by I!i 2 . 14 ^ or ^2.15\} if Ak is a fixed point in O or on dfl, and cr{Ak) is 
defined by or ^2.11 ) if Ak is a limit cycle in O or on dQ. Moreover, if the maximum in 

\2.18\) is attained at exactly one component M := Ak^, then the scaled logarithmic transform 
tCp = —e log Up of the first eigenfunction (normalized by max = 1) converges uniformly in 
O to the maximal viscosity solution W of H2.5\) - /[KU) vanishing on M, i.e. W{x) = dH{x,M). 


3 Passing to the limit by vanishing viscosity techniques 

In this section we pass to the limit, as £ —)■ 0, in equation fl2.2l) and boundary condition fl2.3p 
to get problem fl2.5p - fl2.6l) . We use the standard technique based on the maximum principle 
and the a priori uniform W^’°° bound for obtained by Berstain’s method p, m- 
First, considering fl2.2p at the maximum and minimum points of W^{x) we easily get 

Lemma 2. The first eigenvalue Ap satisfies the estimates minc((r) < Ap < maxc(x). 

Next we establish the bound for ITp in 

Lemma 3. Let Up be normalized by maxup = 1 ( i.e. minITp = 0). Then ||IPe|lwi’°°(o) < C 
with a constant C independent of e. 
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Proof. Following [13] observe that the boundary condition = 0 yields the pointwise bound 


A 

du 


< C'lVhF, 


on dil. 


Therefore, for an appropriate positive function 0 G 




£\ I — 


< 






on dfl. 


(3,1) 


Next we use Bernstain’s method to obtain a uniform bound for u)s{x) = (f){x)\VWe{x)\^, follow¬ 
ing closely the line of |1], Lemma 1.2. In view of fl3.ip either |VB4| = 0 and we have nothing 
to prove, or maxtCg is attained at a point G fl. In the latter case we have = 0 and 


dxjdx 




0|VlF,r <0 at a: = e 


Expanding the left hand side of this inequality we get 
^ , d'^W, d'^We ^ ^ dW, 

2,£(j)a,ij— ^ ^ ^ 2,£cpa,ij 


dxidxjdxk dxk 


^ 4^0,• 


d(j) dW^ d'^cf) 

SOj 


dxj dxidxk dxk 


dxidxj 


|VB4|2. 

(3.2) 


dxidxk dxjdxk 


Using fl2.2p we obtain 

^ dWe ^ d'^We dWe , 3/2 , , 1/2 , . lool 

o o o—- o o o o— + ^{^e +UJ6 + (^e +1 , at x = ^, (3.3) 

OXiOXjOXj^ UX}^ OXk UXiOXj UXf^ \ J 

where we have also exploited the fact that = 0. Substitute now 03.31) into 03.2p to derive 


dxidxk dxidxk 

On the other hand it follows from 02.21) that 




x = f. 


Up < 




d^Wp 


+ 1 


(3.4) 


(3.5) 


dxidxj 

Combining 03. 4 p and 03.5p we obtain Up < C and the required uniform bound follows. □ 

With a priori bounds from Lemma [2] and Lemma [3] it is quite standard to pass to the 
limit in 02.2p . Indeed, up to extracting a subsequence, Wp —)■ W uniformly in O and Ag A. 
Consider a test function $ G C'°°(0) and assume that fU — $ attains strict maximum at a 
point f. Then Wp — $ attains local maximum at fp such that ^ as e —)■ 0. If G O then 
VW,(e,) = V<F(e,) and 


d^Wp 

dxjdx. 


iQ < a. 


^2$ 
dxjdx. 


He) 


and Lemma [2] we get 


and < 0 if G dVL. Passing to the limit as e ^ 0 and using 

f (9$ 1 

/9(V4>(O,O<0 ifeefl, and min ji7(V<F(0,0,^(0 [<0 if ^ e 912. 

Arguing similarly in the case when is a strict minimum of PF — <|) we conclude that IF is a 
viscosity solution of 02.5p - (l2.6p . 
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4 Matching lower and upper bounds for eigenvalues and 


selection of the solution of (l2.5lj —( 12761 ) 


Due to the results of Section [3] we can assume, passing to a subsequence if necessary, that 
eigenvalues converge to a hnite limit A and functions converge uniformly in D to a 
solution W of problem fl2.5l) - fl2.6p as e —)■ 0. In the following four steps we prove that A and 
W(x) are described by Theorem [TJ 

Step I: Significant component(s) of Ah- Recall the dehnition of the partial order relation ^ on 
Ah introduced in [16] as follows 


A!fi,A!' ^ W{A!') = dH{A!',A!) + W{A!). 


(4.1) 


Since the distance function dH{x,y) satishes the triangle inequality and dH{A",A') + 
dn {A', A") > 0 for different components A', A" of the Aubry set Ah, (14.Ih indeed dehnes 
the partial order relation. 

Condition fl 2 . 1 ip assumes that there are hnitely many different components of the Aubry set. 
It follows that there exists at least one minimal component Ai := Ako (such that, WAk 7 ^ A4, 
either M. Ak oi M. and Ak are not comparable). 

Now show that 


W{x) = dH{x,M.) + W{M) in f/ n D, where f/ is a neighborhood of M. (4.2) 


Indeed, otherwise there is a sequence Xi M and a component Ak Ai such that W{xi) = 
dnixi^Ak) + W{Ak)- Then taking the limit we derive W{Ai) = dni-M, Ak) + W{Ak), that is 
Ak A Ai which contradicts the minimality of Ai. 

In what follows a component Ai such that (14. 2 p is satished will be called a signihcant com¬ 
ponent. We have shown that under condition fl2.11l) there is at least one signihcant component 
in the Aubry set Ah- 

Step II: Upper bound for eigenvalues. The crucial technical result in the proof of Theorem [T] is 
the following Lemma whose proof is presented in subsequent four Sections dealing separately 
with four possible cases of the structure of Ai. 

Lemma 4. Let Ai be a significant component of the Aubry set Ah satisfying either i2.1A) or 
1)2.1 A) . Then for sufficiently small 5 > 0 there are continuous functions W^{x), W^fix) and 
neighborhoods Us of Ai such that 


W^^iAi) = 0 and Ws-(x) < W{x) - W{A4) < Wfi{x) m Us nn\Ai, 


(4.3) 


G C‘^{Us n D), uniformly inUsUUL as £ 0 , and 


/ d'^Wt 1 

+ c(&)) > o(M). 


(4,4) 























lim sup 

5^-0 £^- 0 ,^e 


lim sup ( 
■M 


d^Wr 1 N 

(&).5.) + c{i,)) < a{M). 


(4.6) 


Moreover, ifUsHd^l 7 ^ 0 then the functions also satisfy > 0 on UsHd^l, and < 0 
on Us n dQ. 


Now, assuming that we know a minimal component Mi of the Aubry set Ah, we can identify 
the limit A of eigenvalues A^. Consider the difference — Wf~^, where are test functions 
described in Lemma 01 By fl4.3p the function W — Wf — W{A4.) vanishes on M. while it is 
strictly positive in a punctured neighborhood of M.. Then, since converge uniformly 

to IT — Wf as e —)■ 0 in a neighborhood of A4, there exists a sequence of local minima of 
We — such that AT. Moreover, if AT fl Shi 7 ^ 0 then = 0 on (9^2 (locally 

near AT) and therefore G for sufficiently small £. For such e we have 


VVF, = VIT 5 - 


and 


“ij 


d^We 

dxidxj 


— 


d^W, 


(5,6 


dxidx-i 


at X = fe- 


Therefore, 

1 

A, = + -H{VWe{QAe) + c{^e) 

d^Wr 1 

Thus we can use fl4.5p here to pass hrst to the lim sup as e —)■ 0 and then as 5 —)■ 0, this yields 
limsupg^O'^e — o'(AT). Similarly one obtains the matching upper bound so that 


lim Ae = (t(AT). (4.6) 

£->■0 

However, since at this point AT is unknown (it depends on W and thus on the particular choice 
of a subsequence made in the beginning of the Section) equality fl4.6p guarantees only the upper 
bound 

lim sup Ae < max<^ (T(Alfc); Mfc C Ah k (4.7) 

£-i -0 1 J 

where the limsup^^o is taken over the whole family {A^, e > 0}. 

Step III: Lower hound for eigenvalues. Consider a component A of the Aubry set Ah such that 
cr(Al) = max{a{Ak); Ak C Ah}- Introduce a smooth function p{x) such that 


p{x) > 0 in n, p(x) = 0 in a neighborhood of A, and p{x) > 0, when x G Ah \ A 


and consider the hrst eigenvalue A^ of an auxiliary eigenvalue problem 


Q‘^'11^ du^ /I _ 

gaq(a() g^ + (c{x) “ -p(^) )“£ = in 


(4.8) 
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with the Neumann condition ^ = 0 on dil. By the Krein-Rutman theorem the eigenvalue 
is real and of multiplicity one, and Ue being normalized by max^ = 1 satishes Ug > 0 in fl. 
Note that the adjoint problem also has a sign preserving eigenfunction. Then it follows that 

A. < A,. (4.9) 


Indeed, otherwise we have 

Qi^n Qn / 1 \ _ /_ 1 \ 

[c{x)--p{x)jUe-\eUe = - (^Ae - A^ +-p(x) j <0 in hi. (4.10) 

On the other hand, by Fredholm’s theorem the right hand side in (Id.lOh must be orthogonal 
(in L^(hl)) to any eigenfunction of the problem adjoint to fl4.8p . Since the latter problem has a 
sign preserving eigenfunction we arrive at a contradiction which proves fl4.9p . 

Let We '■= —s\ogUe be the scaled logarithmic transform oiue, i.e. Ue = Following 

the line of Section |3] one can show that, up to extracting a subsequence, functions We converge 
(uniformly in hi) to a viscosity solution hFof the problem 


H(VW(x),x) - p(x) =AmQ 


(4.11) 


with the boundary condition ^ = 0, where A = lime_j.o£A£. Note that the argument in 
Lemma |2] yields now bounds of the form — j < A^ < C with some C > 0 independent of e. 
Nevertheless these bounds are sufficient to derive problem (14.lip with the Neumann boundary 
condition. Moreover, since p = 0 in a neighborhood of A one can show that A = 0 using testing 
curves rj from fl2.10p in the variational representation for the additive eigenvalue A (see [7]), 

. 1 rT 


A = — ^im / {L{—v,ri) + p{ri)) df; p solves fl2.8p with 77 ( 0 ) = a; G hlj. 

V 0 

This implies, in particular, that 


W{x) = dniXjA) in a neighborhood of A, 


where dnix, y) is the distance function given by (I2.9p . Then arguing as in second step we obtain 

Ae —!■ (7{A). 

Thanks to (I4.9p this yields the lower bound liminf A^ > max{cr(Alfc); Ak G Ah} complementary 
to fl4.7l) . Thus formula fl2.18p is proved. 

Step IV: Selection of the solution of H2.5\) - I[KU) . Let us assume now that the maximum in 
fl2.18l) is attained at exactly one component A4. Then comparing fl2.18p with fl4.6p we see that 
A4 is the unique signihcant component in Ah, therefore it is the only minimal component of 
Ah with respect to the order relation -<. Thus Ai is the least component of Ah- It follows 
that W{Ak) — W{M.) = dniAk, -M) for every Ak C Ah- Then by fl2.7p the representation 
W{x) = dn^x, M) + W{M) holds in fl. Finally, since mm.-^W{x) = lim^^o min^y IFe(x) = 0 we 
have W{M.) = 0, i.e. W{x) = dnix^M)- 

Theorem [T] is proved. □ 
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5 Construction of test functions: case of fixed points in 


The central part in the proof of Theorem [T] is the construction of test functions satisfying the 
conditions of Lemma H] for different types of components of the Aubry set Ah- Consider hrst 
the case when a hxed point G of the ODE x = b{x) is a signihcant component oi Ah- We 
can assume that W{^) = 0, subtracting an appropriate constant if necessary. Then W{x) is 
given by 

W{x) = dnixA) in a neighborhood U{^) of (5.1) 


We begin by studying the local behavior of W(x) near Consider for sufficiently small 2 ; the 
ansatz 


W{z + ii) = TijZiZj +o{\z\^) 


(5.2) 


with a symmetric N x N matrix T. After substituting fl5.2p into fl2.5p we are led to the Riccati 
matrix equation 

dTQT -TB- B*T = 0, (5.3) 


where Q = [aij{^) 






i,j=l,N 


Next we show that fl5.2l) holds with T being the maximal symmetric solution of fl5.3p : for 
existence of such a solution see, e.g., [12] or [1] . To this end consider the solution D of the 
Lyapunov matrix equation 


D{4:TQ -B) + (dTQ - R)*D = 21 


(5.4) 


given by 

D = 2 f (5.5) 

J —00 

By Theorem 9.1.3 in [12] all the eigenvalues of the matrix dQT — B have positive real parts, so 
that the integral in (15.5p does converge. Set 

Tf = T± 6D. 


Then T^ satishes 


dT.-gr,- - TJB - B*TJ < -61 


(5.6) 


for sufficiently small 5 > 0. 

Introduce the quadratic function Wf (x) := Tj{x — ^) ■ {x — . Thanks to fl5.6p this function 

satishes 


H{VWi{x),x)<-^-\x-^\^ 


in a neighborhood of 


(5.7) 


This yields the following result whose proof is identical to the proof of Lemma 16 in [16] (see 
also the arguments in the proof of Lemma 0 below). 
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Lemma 5. The strict pointwise inequality {x) < W{x) holds in a punctured neighborhood 
of ^ for sufficiently small 5 > 0. 

Next consider the function Wf'^x) := T'g {x — ■ {x — 

Lemma 6. The strict pointwise inequality (x) > W(x) holds in a punctured neighborhood 
of ^ for sufficiently small h > 0 . 

Proof. According to fIS.ip . the following inequality holds 

W{x)< f L{-v{t),^ +7](t)) dr 

Jo 

for every control v{t) such that the solution of the ODE 

i){t) = v{t), ?7(0) = z := X — 


vanishes at the hnal time t and remains in a small neighborhood of 0 for any 0 < r < f. We 
can take the hnal time t = +oo and construct n(r) by setting n(r) = —{dQT — B)ri{T), where 
r]{T) in the solution of the ODE 

T] = -(4Qr - B)r], r]{0) = z. 


As already mentioned (see Theorem 9.1.3. in lEl) all the eigenvalues of the matrix AQT — B 
have positive real parts, therefore \ri{T)\ < C\z\ and pir) —)■ 0 as r ^ +oo. Moreover, the 
latter convergence is exponentially fast. 

Thus we have 

L(-n(r),r/(r)+0 = if + + Hh)){-Vi + Hv)) = 

= + BikPk){-Pj + Bjipi) + 0{\p\^), 

where a^^{x) denote the entries of the matrix inverse to (aij{x)). Next recall that Qij = 

aij{^) and that T solves (15.31) . Taking this into account we obtain 

noo noo 

/ L{h + i, -v{r)) = - + BikVk)i-Vj + Bjipf) + 0{\z\^) = 


= -2 


Tp-r]dT+ / Tp ■ {p P Bp)dT + 0{\z\^) = 


= Tz- z + 


p ■ (-4rgr + te + B*T)pdT + = 


= TijZiZj + Od^d). 


Finally, since by the dehnition of T^, T^ = T + 6D with D > 0 , then for sufficiently small z ^ 0 
we have 

W{z + 0 < Tz • ^ + Od^n < Tjz ■ z. 

□ 
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Lemmas O and [H] show that functions do satisfy conditions of Lemma 01 To complete 
the proof of Lemma |4] in the case of M. being a hxed point in we dehne functions simply 
by setting := W^. Thanks to fl5.7p we have 




/ u i \ 

hrnsup(a,,(e.)^^(e.) + + c(e.)) < -2a,,(0(rp- - <5A,) + c(0, (5.8) 


as soon as when e —)■ 0. According to Proposition 20 in [16], —2aij{^)Tij + c{^) = cr({.^}), 

thus fIS.Sp yields fld.Sp . Similarly one verihes that satishes fl4.4p . 


6 Construction of test functions: case of fixed points on 

dn 

Consider now the case of signihcant component of the Aubry set Ah being a hyperbolic hxed 
point ^ of the ODE x = br{x) on dQ, where br{x) denotes the tangential component of the 
vector held b{x) on dfl. As above, without loss of generality, we assume that IT(^) = 0 . 

It is convenient to introduce local coordinates near dfl so that x = X{zi,..., z^) with 
zn = being the distance from x to dfl {zh{x) > 0 if x G D) and z' = {zi, zh-i) 

representing coordinates on dfl in a neighborhood of the point The latter coordinates are 

chosen so that the map X{z\ zn) is C^-smooth and z'{^) = 0. Moreover, the matrix _ 

is orthogonal when z' = Q and z^ = 0 (at the point In these new variables equations fl2.5p 
and fl2.2p read 

S{V,W,z) = 0 (6.1) 

and 

- ea^(X{z))%-'(z)X(r-i(,)?^X S(v,w„z) = e(X,- c(X(z))), (6.2) 

where 

S{p,z) = aij{X{z))%~i\z)ri~\z)pkPi - bi{X{z))rAi\z)pk 

and ( TX^iz )) _is the inverse matrix to ( ^^(z )) _. Note that according to hypothesis 

V / i,j=l,N J i,j=l,N 

bi{X{z))Txl{z) < 0 for sufficiently small \z\. (6.3) 

Like in Section |5] we construct the leading term of the asymptotic expansion of W near the 
hxed point ^ in the form of a quadratic function. Taking into account the boundary condition 
= 0 (that is 02.61) rewritten in aforementioned local coordinates) we write down the following 
ansatz 

W{X{z\zn)) = fijz'iz'j + o{\z\^ + zIj). 

with a symmetric [N — 1) x [N — 1) matrix T satisfying the Riccati equation 

4fQf -TB- B*f = 0, (6.4) 
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where Q = [aij{O%i\0)Tij ^(0) 


k,l=l,N-l 


imdB= r„'(0)|^({)S(0) 


Note that 


k,l=l,N-l 


B is nothing but the matrix in the ODE z' = Bz' obtained by linearizing the ODE x = br{x) 
near ^ in the local coordinates z' = {z [,..., z'^_i) on dVL. 

Let r be the maximal symmetric solution of fl6.4p . and let D be a solution of the Lyapunov 
matrix equation 

b{AfQ-B) + {AfQ-Byb = 2I. (6.5) 

By Theorem 9.1.3 in [12] all the eigenvalues of the matrix dLQ — B have positive real parts, 
therefore fl6.5p has the unique solution b given by 

^0 


D = 2 


g(4rQ-S)*tg(4rQ-S)t^^^ 


which is a symmetric positive dehnite matrix. Now introduce functions 


Wy{z', zn) = (f ± 6b)ijz'z' ± 6z 


'N 


( 6 . 6 ) 


depending on the parameter 5 > 0 . 


Lemma 7. Let 6 > 0 be sujficiently small. Then, for small \z\ y 0 such that X(z) G D, we 
have 


Wi(z) < W{X{z)) < Wf(z). 


(6.7) 


Proof. By virtue of the dehnition of it suffices to prove (16.711 with non strict inequalities 
in place of strict ones an then pass to slightly bigger 6. 

The proof of the inequality Wf < hh is based on the following two facts. First, we use the 
fact that W{x) = dnix,^) in a neighborhood of f. Moreover, for a given > 0 there exists 
(5 > 0 such that if |a: — ^| < 6 then the minimization in (12.9p is actually restricted to testing 
curves rjb) which do not leave the set {|?7 —■Cl < (otherwise arguing as in [T^ Lemma 19] one 
can show that f is not an isolated point of the Aubry set Ah, contradicting (12.1111 ). Second, 
considering, with a little abuse of notation, Wy{x) = Wy{X~^{x)) we have for sufficiently 
small (5 > 0 

o dwr 

H(VWg , x) < —6\x — in D, and ^ = 0 on dQ, (6.8) 

when |x — .^1 < 6' with some 5' > 0 independent of 6. This follows from the construction (16.6p 
of Wy and (16.4p . (I6.5p . also taking into account (16.3p . 

Assume that |x — ^| < 5, and let r]{T) be a solution of (12.8p satisfying 77 ( 0 ) = x, rjit) = f 
with a control v{t) such that \ri{T) — .^| < 5' for all 0 < r < t. Then 

VWs-{v)-vdT= [ VWs-{v)-{-v{T))dT, 

Jo 



where we have used the fact that = 0 on dVt. It follows by Fenchel’s inequality p ■ {—v) < 
L{—v, 77 ) + H{p, 77 ) that 


Wf{x)< / L(- 7 ;, 77 )dr+ / H{XWs-,p)dT< / L{-v,p)dT. 
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Therefore by fl2.9l] we obtain Wf{x) < W{x). 

In order to prove the second inequality in fl6.7p for a given x = X{z',zn) we construct a 
test curve p(r) hrst on a small interval (0, At) by setting rj^r) = A(z', CAr(r)), Civ('r) being the 
solution of ODE Cn{t) = bi{X{z'XN))Tj^i{z'XN) with the initial condition Civ(O) = zn, and 
choosing At from the conditions Cv(At) = 0, (nX) > 0 for r < At. Thanks to fl6.3p we have 
At = 0{zn)- Then, since 

m = ^X'XN)bMTxXz'XN) 

OZn 

BX f)X 

= -^X'XN)bXXr-,\z'X n) - -^{Z'XXUXT-XX'XX = hX) + 0 {\z\) 

^ 3 

(recall that the tangential component b^- on dQ vanishes at the point we obtain 

j-At 


(6,9) 


Next we construct riX) for r > At which connects point X{z',0) to Following closely the 
line of Lemma [6] we introduce CX) by solving the equation C'X) = —{AQT — B)(' with the 
initial condition = z' and set rjX) = ^X'X), 0). Then rjX) solves (12.81) for r > At with 

dX 

vX) := XX + XvXuX) - X'T'mXXijiX'TijXX'^imCm 

(note that ^(C',0) = and b^X) + > 0 as soon as \z'\ is 

sufficiently small) and using (16.4p we obtain 


'At 


L(-n(r),r 7 (r))dr = - / a*^(0(-'y* + dr + 0 (| 2 ;|^) 

4 J At 


= 4 


'At 


/ rlX r)X 

a^Xi){-^^X)Qutu'm + ^X)rxlXXXi)rJX)^mnC 

/f)X ~ f)X ~ \ 

X [-^X)QklflmCm + -^X)rXXXaklX)rJiOXmnCn) + 0(|z| 


= 4 
= -2 


At 


fc'-gfc'dr + 0(|z|3) 


'At 


fC',C'dT+ / fC',(C' + BC')dT + 0(|s|=>)=f„,4j'+0(|z|=>). 


'At 


The required upper bound W < now follows from (16.9p and (16.101) . 

Thus functions satisfy conditions of Lemma HJ moreover it follows from 
junction with (16.4p . (16.5p . taking also into account (16.3p . that 

S(VzW^X)yZ) A 0 and S(V zWf X)y z) < 0 when \z\ is sufficiently small. 

Then we set 

iy± (/, zn) = W^X', zn) T sXn, 

and verify (similarly to the case of interior hxed points) that conditions (14. 4 p and 
satished. 


( 6 . 10 ) 

□ 

in con- 


are 
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7 Construction of test functions: case of limit cycles in 

Q 

We proceed with the case when a signihcant component of the Anbry set Ah is a limit cycle, 
assnming hrst that it is sitnated entirely inside Namely, let ^{t) be a periodic solntion of 
the ODE ^ = b{^) whose minimal period is P > 0. We assnme that C = {^{t) \ t E [0, P)} C D, 
h{x) 7 ^ 0 on C and C is a hyperbolic limit cycle, i.e. the linearized Poincare map associated 
to this cycle has no eigenvalnes on the nnit circle. In order to stndy the local behavior of W 
near the cycle C, perform a O^-smooth change of coordinates x = X{zi, ..., zn-i, zh) with z^ 
representing the arc length along the cycle and z' = (zi,..., zn-i) being some hxed Cartesian 
coordinates in the hyperplanes orthogonal to the cycle. Also we assnme that C is oriented by 
the tangent vector ^(0/1 ^(01) P = 0 on C. With this change of coordinates eqnations 
fl2.5p and (12.21) take the form similar to (16.11) and (16.2p . Assnming as above that W{C) = 0, we 
postnlate in the vicinity of the cycle (for snfficiently small \z'\) the following ansatz for W: 

W{X{z\zH))=T,,it)z:z'^+o{\zf), (7.1) 

where t refers to the parametrization of the cycle determined by the eqnation ^ = b{^), t G [0, P). 
Snbstitnte W in (16.11) to hnd, after collecting qnadratic terms and neglecting higher order terms, 

r = 4rgr-pp- pT, ( 7 . 2 ) 

where Q{t) and B{t) are P-periodic {N — 1) x (A^ — 1) matrices whose entries are given by 

%,(() =««KW)7;7'(0,zkK(()))7;-'(0,zjvK(())). (7.3) 

_ au av r\ / BX \ 

Bu{t) = (0, (^(0) (0, ^„K(()))-ry ■ (0, (^(o, ^a(«())) j. 

(7.4) 

Recall that TiJ^ [z', zn) denote the entries of the matrix inverse to {^^{z', zn)) ^ 
brevity abusing slightly the notation we set 

r)X 

Y\t)-.= TZ(Q,z^(m), 7I,(() = .^(0,^aK(())). (7.6) 

The matrix Q{t) being positive dehnite, it is known P that Riccati equation (17.2p has 
a maximal symmetric P-periodic solution P(f). We next show that (17. ip does hold with the 
mentioned maximal solution P(t) under our standing hyperbolicity assumption on C. Note that 
the ODE z' = Bz' corresponds to the linearization of i = b{x) on the cycle C written in local 
coordinates; thus assuming the hyperbolicity of C we require that the fundamental solution of 
the ODE _ 

^(f,r) = P(t)<l>(f,r), $(r,r) = /, 

evaluated at t = r -|- P, has no eigenvalues with absolute value equal to 1. 
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Lemma 8. The following bound holds uniformly in t E [0, P) for sufficiently small \z'\, 

w{x(z',ZM(m))<’ri,(tW, + c\z'\Hogf. 


(7.6) 


Proof. We make use of variational representation fl2.9p . A natural guess about optimal test 
curve in fl2.9p is that its first — 1 local coordinates are given by (cf. Sections [5] and [6]) 

C'(r) = (B- 4g r)C'(r), r > t, C'H) = z'. (7.7) 

Thanks to Theorem 5.4.15 in [1] the solutions of fl7.7p are exponentially stable, i. e. |CX'^)I — 
Ce~^'^\z'\ for some 5 > 0. The choice of the last local coordinate is a little bit involved. We set 
/^(r) := A(g(r), 2 Ar(.^(r)) + (CAr(r)) and want to choose CAf(T) in such a way that |Cv('r)| < C\z'\, 
and 

hifr)) - =4aij(^(r))7;"^(0,^(r))r;^(r)C(r) + 0(|/H. (7.8) 

We skip for a moment the proof of the existence of such C,n- It will be given later on. 
Considering fl7.8p we obtain 


L(-r),r7)dr = - / a*^(r7(r))(-?7i(r) + 6i(r7(r)))(-r)j(r) + g(r7(r)))dr 


< 4 / (r(r)C'(r)) • (g(r)r(r)C'(r))dr + CT\z^ 


/|3 


In view of (ThTP and (TT^ we have 


L{-f r/)dr < / (TCO ■ {PC' - Odr + CT\z 


/|3 


dr 


(rg ■ c')dr + / (rg ■ {bc,' + c') + rg ■ C')dr + CTf 


/|3 


< T{t)z' -z' + {T - dV QT + T B + B* r)C' ■ C')dr + CTf' 


,/|3 


= T{t)z' ■ z' + CTf 


/|3 


If we choose T := -f log-r^, then dist(r 7 (T), C) = 0{\z'f), and 


L{—r], ri)dT < T{t)z' ■ z' + C\z'\'^ log — 


(7.9) 


For constructing (n we will need the following facts. From the dehnition of fj in fl7.5p it 


follows that Tmir) = g(r)/|^(r)|. Then, since g(r) = £^(^(r))^j(r), we have 


*■ iUrWMr) - TMr) = = 7-«(r)Mh. « = 1..... ]V. 


dx 


le(r)|^ 


le(r)P 
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Multiplying this by T ^(r), we conclude that 




iffc = iV, 

0 iik<N. 


We proceed with constructing Cat. By the dehnition of rj we have 


(7.10) 


Vi{r) - biivir)) = ^ (^^*(r) + 7'fc(r)CUr) + 7Iiv(r)C7v(T)) - Hvir)) + 0{\z'\^) 

= (fik{T) - ^UM)7jk(r)^(k(^) + ('^v(r) - ^(^(r))7'Ar(r))Cv(r) 

+ 7I.(r)C^(r) + TMrKN(T) + 0{\zf) 

= Tik{T)%~\T)i^(frk{T) - ^(^(r))7^fc(r))Cfe(r) + ('^v(r) - ^(^(r))7^v(r))Cv(r)} 
+ 7ifc(r)C(.(r) + 7iv('7‘)Cw('7‘) + 0(|2;'|^). 

(7.11) 

Substituting for (^' the expression on the right-hand side of fl7.7p and considering fl7.10p yields 
Hvir)) - Vi{^) = 4aij(^(r))7;7^(0, ^(r))ri^(r)C™(r) + 0{\zf) - 7l7v(r)C7v(r) + TiN{T)R{T), 
where 


Rir) = 


^(^W) (^7^v(r)Civ(r) + 7'i(i‘)C)(7-)^ 

(TiN{T)CNir) + rii{T)Ci{T)^ -^T^lir)aij{^{T))rij\T)TirniT)Cir) 


Thus, in order to make fl7.8l) hold, we choose Cw('r) as a solution of the following equation 

Cw(r) = i?(r) (7.12) 


with the initial condition Cv(^) = 0. 

From fl7.1Up it follows that |^(r)| solves 

4|?(r)| = (r„-'(r)^(e(r))7}A.(r) - r„-/(r)7-A-(r))|e(r)|, 

and we can write the solution Cn{t) of (I7.12p as 

T 

Cn(t) = |«t)| j (^r,^.i(i,)^K(ii))7-,(ii)C,'(s) - t„1(s)%(s)Cm'^ ^ 

“ , ‘ (7.13) 

-4|«t)| / r„-‘(s)ayK(s)) 7 ;-‘(s)r,„(s)C(ii)^. 

J l^(s)| 

From (I7.13P we derive the uniform bound |CAr('r)| <C\z'\. 
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It remains to construct ri{-) for r > T in such a way that it reaches the cycle in a hnite 
time. To this end we set 

,,(t) = x{('(T)(T +1 - r),+ u(r)l«(T)l/l«(r)l). r < T < T +1. (7.i4) 

Then for every i = 1,... ,N 

rW) = T({.(r) + 7-„(r)&(r)|«r)|/|«T)|) + 0{\zf) 

= Ur)+Ur)UT)/K(T)\+0(\zf) 
f)h ■ 

= hiar)) + ^(e(r))0(r)Cv(T)/|e(T)| + 0(|zf) = k{v{r)) + 0(|zf). 
Therefore, the following bound holds 

T+l 

J a'^{v{r)){-Vi + bi{v))i-Vj + bj{v))dT < C|z'|h (7.15) 

T 

Combining the last relation with fl7.9p yields fl7.6p . □ 

In order to construct a sub- and supersolutions of fl6.ip we consider the solution D{t) of the 
matrix equation 

D+^{B-AQT) + {B-AQVyD =-21, (7.16) 

given by 

OO 

D{t) =2 J A/*{T,t)A/{T,t)dT, (7.17) 

t 

where T (r, t) is the fundamental matrix solution of 

= (B(t) - 4Q(T)r(r))^, ^((,0 = /. 

As already mentioned in the proof of Lemma [H] this solution T(r, t) decays exponentially as 
r -l-cxo and therefore the integral in fl7.17p converges. Then it dehnes a P-periodic positive 
symmetric solution of fl7.16p . It follows from fl7.2l) and fl7.16p that T^ := T ± 5D satisfy for 
sufficiently small <5 > 0 

4rt or; - dr; - r; b - bt; > si and 4r; or; - dr; - r; b - bt; < -si. 

Now dehne the functions Wy{z) by 

where zn = ^ 7 v(^(t)), 

these functions satisfy 

A A 

S{'VzWy{z),z) <--Iz'l"^ and S{VzWy{z), z) >for sufficiently small |P|. (7.18) 
The latter inequalities follow directly form the dehnitions of Wf{z) and W^^z). 
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Lemma 9. For sujficiently small 5 > 0 the strict pointwise inequalities (z) < W{X{z)) < 
W^{z) hold for z' from a punctured neighborhood of zero. 

Proof. The first inequality Wf~{z) < W{X{z)) can be proved similarly to Lemma 19 in [16] (see 
also the proof of Lemma [7|), using fl7.18p . The second inequality W{X{z)) < W^{z) follows 
immediately from Lemma [HI □ 

At this point we have constructed functions satisfying conditions of Lemma jH Next 
we dehne the test functions by 

Wt - Wf - 

where z^ and t are related by z^ = ZN{^{t)), and ^f{t) are periodic solutions of the ODEs 

p p 

(t) = -2tr(Q(t)r^ (f)) + c(^(f)) + p y tr(g(r)rf (r))dr j c(^(r))dr. (7.19) 

0 0 

The hrst two terms on the right-hand side here are introduced in order to compensate the 
discrepancy of order e in equation fl6.2p . Indeed, the test functions constructed in this way 
satisfy for sufficiently small \z'\ 

d ( dW^\ 1 

±aij (X(z)) ^ki(z)-^(^aij(z)-^J T -S(V,W^% z) T c(X(z)) 

p p 

<yy c(^(r))dr-y j tr(g(r)r^(r))dr + 0(£+Iz'l). 

0 0 

In order to complete the proof of the fact that satisfy fl4.4p and 04.51) it remains to observe 
that tr(g(r)rf (r))dr ^ tr(g(r)r(r))dr as 5 —?■ 0 and use the identity 

2 rtr(g(r)r(r))dr = E log 0i 

ei>i 

(see Proposition 5.1 in im). where 0j are absolute values of eigenvalues of the linearized 
Poincare map (corresponding to the ODE i. = b{x) near C). 

8 Construction of test functions: case of limit cycles on 

dn 


In the case when ODE x = br{x) on dQ has a limit cycle C which is signihcant component 
of the Aubry set, the analysis combines the ideas of Section [6] and Section [71 We pass to the 
local coordinates in a neighborhood of C via a map x = X{zi ,..., zn-i, zn), where zn = zn{x) 
is the distance from x to dVL (positive for x E fl) and {zi,... ,zn-i) are coordinates on dVL. 
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The coordinate Z]^-i{x) represents the arc length parametrization on C and other coordinates 
z' = {zi,, zj^_ 2 ) are chosen so that the map X{z', zjsi-i, zn) is C^-smooth, moreover z' = 0 
when X E C, and ^ orthogonal matrix when Zn = 0 and z' = 0 (on the 

cycle). This change of coordinates leads to eqnations of the form (16.Ij) and (16.2^ for W{X{z)) 
and W,{X{z)). 

We use the following ansatz for W, 

W(X(z))^Y(t)z:z'i+o(\zf), ( 8 . 1 ) 

where T is now {N — 2) x {N — 2) symmetric P-periodic matrix (P being the period of the 
cycle C), and t refers to the parametrization t ^{t) of C such that ^{t) = bT-{^{t)). Moreover, 
T is chosen to be the maximal P-periodic solution of the Riccati matrix equation 

—f = dfQT - fP - P*f, 
dt 

with (N — 2) X (N — 2) matrices Q{t) and B(t) whose entries are given by the same formulas 
as (Q and dZai). 


Lemma 10. For sufficiently small \z'\ and \zn\ the following hound holds uniformly int E [0, P) 

W{X{z',ZN-i{f{t)),ZN) < fij{t)z-Zj + C{\z'\^ \og^^F \zn\\z'\^ + kArH- ( 8 . 2 ) 

Proof. First consider the case z^ = 0. As in Lemma |H] we use representation 02.91) and consider 
the solution Ci'^) of fho ODE C{t) = {B — 4Qr)(^(r) for r > t with the initial condition 
C'(t) = z' . It decays exponentially as r —?> oo, |C'| < for some 5 > 0. Next we 

introduce Cv-i analogously to Cw introduced in Lemma [H i.e. Cn-i solves 

Cv-i(r) = 7^)v-i)i(0>^^-i(^(^)>0)^(^(^))('7i(iv-i)(r)C7v-i(r) + 7^z(r)C;(r)^ 

- 2;jv-i(^(r)), O) (^7;(jv-i)(r)Civ-i(r) + '^z(r)C/(r)^ 

-47()^^_i)i(0,Z7v-i(^(r)),0)ay(^(r))7;-^(0,z^_i(^(r)),0)R^(r)C(r), r > t, 

where is the matrix inverse to (|J(2;))ij=Tyv and 7Ij(r) = |J(0, . 2 Ar_i(^(r)), 0). 

Finally we dehne "gfr) by 

^v-i(^(r)) + Cv-i(r), 0), t<T <T 

^ [X(C'(T)(T + 1 - r), zw-i(e(r)) + Cv-i(T)|e(r)|/|e(T)|, 0), T < r < T + 1, 

with T := I log -j^, and the control n(r) by 

fr) +i/(r/) ( 6 ^( 77 )+Pi(r)), t<T<T 
v{t) = < 

\i) + n{r])b^{r]), T<t<T+1, 
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where zn_i = and 

^l(^) = -7^^^_i)i(0,^7V-l,0) (fi^N-l){T)CN-l{T)+filiT)CiiT)+Aaij{^{T))Ti~^{0, ZN-uO)f i^C'm) ■ 
Letting 

, . i Mv) + Ri{r), t<T <T 
a(r) := < 

\^K{r]), T < T <T + 1, 

observe that for this control v{t) the pair (//(r), a(r)) solves fl2.8p on (t, T + 1) with the initial 
value r]{t) = x{= X{z',^{t),0)), as far as a{T) > 0 for all r G {t,T + 1). Since > 0, 

the latter condition is satished, provided that \z'\ is sufficiently small. Then the proof of fl8.2p 
follows exactly the line of Lemma El 

In the case when zn{x) > 0 we construct a curve rj^r) connecting x with a point y on dVt 
by setting 


7]{t) = X{z', ZN-i{i{t)), zn + — r)) for all r > f such that zn + — r) > 0. 

Let f + At be the time when r]{T) reaches dQ (at the point y = r]{t + At)) then At = 0{zn)- 
It follows from the construction of ? 7 (r) that 

t+At 

{'n{T)){-r]i + bi{T)){-f]j + bj{r]))dT < C{\z'\‘^ + z%)zn- 
Then extending ri{T) along dQ as described above we complete the proof of the Lemma. □ 


Now we construct test functions W^{z', Z]\f-i{^{t)), zn) := (T ± 6D)ij{t)zlz'j ±6zjf for (suf- 
hciently small) <5 > 0, where the P-periodic symmetric matrix D{t) >0 is dehned analogously 
to fl7.17p . Then 

^(V.IT^-,^) < -6{\zf + Km)zN) 

for sufficiently small \z'\ and zn > 0. This yields the following bound 


IT,- < W{X{z)) for sufficiently small \z'\ and zj^ (when |P| + > 0) 

whose proof is analogous to that of the lower bound in Lemma [71 Thus functions IT,^ satisfy 
the conditions of Lemma jH Finally, we dehne the test functions by 

■= ± T ^‘^zn, where zn-i = ZN-i{^(t)), 

with being solutions of 


p p 

= -2tr(Q(f)f^(f)) + c(^(f)) + jj tr(Q(r)f^(r))dr j c(^W)dr. 

0 0 

These functions IT,^^ satisfy the conditions of Lemma jH 
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